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Primordial black holes and the observed Galactic 511 keV line 
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The observed 511 keV line from the Galactic Bulge is a real challenge for theoretical astrophysics: 
despite a lot of suggested mechanisms, there is still no convincing explanation and the origin of 
the annihilated positrons remains unknown. Here we discuss the possibility that a population of 
slowly evaporating primordial black holes with the mass around 10^® — 10^^ g ejects (among other 
particles) low-energy positrons into the Galaxy. In addition to positrons, we have also calculated 
the spectrum and number density of photons and neutrinos produced by such black holes and found 
that the photons are potentially observable in the near future, while the neutrino flux is too weak 
and below the terrestrial and extra-terrestrial backgrounds. Depending on their mass distribution, 
such black holes could make a small fraction or the whole cosmological dark matter. 
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Introduction — It is now clear that in the central 
region of the Galaxy electron-positron annihilation pro- 
ceeds at a surprisingly high rate. Confirming precedent 
measurements [s^, the SPI spectrometer on the INTE- 
GRAL satellite has detected an intense 511 keV gamma 
ray line flux (Bulge component) d, [1] 



$511 = 1-07 ± 0.03 • 10"^ photons cni"^ s" 



(1) 



with a width of about 3 keV, consistent with two di- 
mensional gaussian distribution aligned with the Galactic 
Center and with a full width at half maximum (FWHM) 
of about 8°. More recently, the SPI spectrometer has 
also provided evidence for the disk or halo component Q 
of this line. 

Non-relativistic positrons in the interstellar medium 
can either directly annihilate with electrons or form 
positronium. In the first case they produce two 511 keV 
photons, while in the second case the annihilation chan- 
nels are different for para-positronium (formed with 25% 
probability) and ortho-positronium (formed with 75% 
probability). Para-positronium decays into two 511 keV 
7S, while ortho-positronium decays into three photons 
with continuous spectrum. From the relative intensi- 
ties of the 511 keV line and the three 7 continuum, one 
can deduce that the total annihilation rate is 3.6 times 
larger than the one arising from consideration of flux ([1]) 
only 0. Thus, assuming that the Solar System lies at a 
distance r = 8.5 kpc from the Galactic Center, we can 
conclude that the total rate is about 3-10'*^ annihilations 
per second. 

From the theoretical point of view, the issue is to 
identify the source of these galactic positrons. Many 
production mechanisms have been suggested, but so far 
no one is completely satisfactory. Among conventional 
mechanisms, which do not demand new physics, there 
are type la supernovae 0], low mass X-ray binary sys- 
tems 01 (see also [1|) and the energetic electrons and 
photons created by accretion on the super-massive black 
hole at the Galactic Center [§]. A similar mechanism of 



positron production by collision of energetic photons pro- 
duced in accretion to super-massive central black hole 
and to surrounding primordial black holes with mass 
about 10^^ g was considered in ref. [li|. More exotic 
scenarios include annihilating light dark matter parti- 
cles [ll|, deca yin g unstable relics and, in particular, ster- 
ile neutrinos [l^l, MeV right-handed neutrino interact- 
ing with baryonic matter |13'|, strangele ts fl4l| . positrons 
originating from primordial antimatter [l5| , from decays 
of milli-charged particles [l6|, and possibly a few more. 
However, type la supernovae have a different galactic dis- 
tribution and their rate is probably an order of magnitude 
smaller than the one necessary to explain the observed 
flux 0,113 (see however (H). Low mass X-ray binaries 
also do not fit the enhanced Bulge component 0] . Lastly, 
light dark matter particles are theoretically not well mo- 
tivated, might be inconsistent with observations 19], and 
their mass should be very close to the electron one, be- 
cause from the comparison of the Galactic gamma ray 
emission above and below 511 keV, one can conclude that 
the injected energy of the positrons cannot be larger than 
about 3 MeV 0]. 

In this letter, we discuss the possibility that positrons 
are produced by evaporating primordial Black Holes 
(BHs). As in the case of MeV dark matter, some fine 
tuning is needed too, namely, the BH mass distribution 
should be peaked in the interval 10^^ — 10"'^^ g. The idea 
that the primary source of low energy positrons in the 
Galaxy could be primordial BHs was first considered in 
ref. [23| . The picture was later discussed in ref. |21| , with 
the conclusion that primordial BHs could unlikely be re- 
sponsible for the 511 keV line from the Galactic Center, 
unless the primordial BHs are more strongly clustered 
in the halo than the other halo material. However, the 
authors of ref. [2l| assumed that the initial mass dis- 
tribution of BHs, dN/dM, is scale invariant, while we 
have considered a mass distribution which is peaked at 
a particular value. In this case primordial BHs can ac- 
count for the observed 511 keV photons and, at the same 
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time, could make even the whole cosmological dark mat- 
ter, with no contradiction with the observed gamma ray 
backgrounds and gravitational lensing data, if the mass 
spectrum of BHs has a pronounced maximum in the in- 
terval 10^^ — 10^^ g. A mechanism of creation of primor- 
dial BH dark matter with a peaked mass distribution was 
suggested e.g. in ref. [2^ . 

Primordial BHs: general features — It is 
well known that at the semiclassical level BHs are no 
longer one-way membranes, but emit the Hawking ra- 
diation [2^. Restricting to the simplest case of a 
Schwarzschild BH with mass AI, the emission rate of 
the particles of species i with the energy in the range 
{E,E + dE), orbital momentum I, third component of 
the orbital momentum m and polarization s is [33 |: 
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Urns 



dE 



dt 



(27r) exp(£:/r) ± 1 ' 



(2) 



where Tums is the so-called graybody factor equal to 
the absorption probability for an incoming wave with the 
specified quantum numbers, T is the Schwarzschild BH 
temperature 
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and the signs ± are for fermions and bosons respectively. 
At high energies {GnME > 1), ^iims oc {ME)^, since 
the cross section for each kind of particle approaches the 
geometrical-optics limit, and BH radiates basically as 
black body of the same temperature. 

As one can easily see from eq. ([3]), the temperature of 
an astrophysical BH with the mass of the order of the 
Solar mass Mq « 2 ■ 10'^^ g, or larger, is very low, below 
10~^ K, and the corresponding particle emission is com- 
pletely negligible. However, lighter BHs could have been 
produced in the very early universe by a lot of possible 
mechanisms (for a review of production mechanisms and 



observational bounds see e.g. [2j| and references therein) 
and today evaporate at a much higher temperatures. 
From eq. ([2]) , one can compute the total particle emission 
rate and the total BH mass loss rate |25i, |26j| 
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= -4.8 • 10 



exp{E/T) ± 1 (27r) G%M^ 



20 



(4.5 • 10-4) 
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MeVs^\ 
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where a' and a are numerical coefficients depending on 
the BH mass M and the particle content of the theory. 
In the Standard Model, for M ^ 10^^ g, one finds a' = 



1.6 • 10-3 and a = 2.8 • 10"^, while, for 5 • lO" 
10" g, a' = 2.6 • 10-^ and a = 4.5 • 10"* ^ 
eq. (O, we can find the BH lifetime 



< M < 
5| . From 
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where = a{Mi) and Mi is the initial BH mass, since 
the BH spends most of its life near its original mass. 
From eq. © we see that all the primordial BHs with the 
original mass lighter than 5 • 10^^ g have already evapo- 
rated. 

Primordial BHs in our Galaxy — Let us start 
with a first estimate, just to show that primordial evap- 
orating BHs can be viable candidates to explain the ob- 
served 511 keV line from the Galactic Bulge. For this 
purpose, we assume that all the BHs have the same mass 
AI — 6 ■ 10^^ g. The positron production rate per BH is 



dN, 



dt 



0.7 



GnM 



2 • 10 



18 e+ s- 



(7) 



where a'^^ = 0.49 ■ 10"'^ and the factor 0.7 is the cor- 
rection due to the finite value of the positron mass. In 
order to explain the observational data, which suggest 
that the positron production rate is 3 • 10*^ e"*" s~i in- 
side a spherical region of the radius r = 500 — 700 pc, 
we need about 1.5 • 10^^ primordial BHs, whose total 
mass is ~ 4.5 • 10® Mq. Here positrons are injected into 
the Galaxy non-relativistically, together with electrons, 
photons and neutrinos (and possible other unknown light 
particles). 

The neutrino flux on the Earth would be extremely 
weak, about 0.02 i' cm~^ s~i, impossible to be detected 
by present and foreseeable future experiments. 

The photon flux from the Galactic Center would be 
about an order of magnitude less intense than that of 
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1: Gamma ray spectra from primordial BHs with mass 
6 • 10^^ g (blue solid curve) and of the diffuse background 
(red dashed curve) in 7 cm~^ s~^ as a function of energy E 
in MeV. The number of BHs is normalized by the condition 
that they produce the observed positron flux. 
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neutrinos [36|, i.e. 

$7,BH « 1-8 • 10"^ 7 cni"2 



(8) 



and may be observable as a possible bump in the diffuse 
background of the Galactic Center. The presence of the 
bump depends on the exact position of the peak in the 
photon spectrum from evaporating BHs. The correct BH 
photon cross section depends also on the energy of the 
emitted photon, that is a-y = (M, i?) , and has to be 
computed numerically [25l . [26j . In this paper we use the 
photon cross section that one can deduce by fitting the 
curve reported in Fig. 1 of ref. [2^ . The gamma spectrum 
from primordial BHs is plotted in Fig. |1] together with the 
measured Galactic continuum one M, u3 
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E 



0.511 MeV 



-1.75 
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where A = 7 • 10 7 cm ^ s ^ MeV ^ is the normaliza- 
tion factor dX E ^ 0.511 MeV. 

We can also notice that the cosmological dark mat- 
ter may be made of 6 • 10^^ g primordial BHs: their 
number density today would be ubh = 4 • 10^'*^ cm^"^ 
and they would have created a cosmological MeV gamma 
background n-y ~ 1.1 ■ 10~^^ cm~'^ or, equivalently, an 
isotropic flux <I>^ « 0.3 7 cm^^ s~^. The latter is about 
what we observe, so it is not unreasonable that the Dark 
Matter (DM) in the Universe is made of BHs with mass 
M = 6-101^ g. 

If the mass fraction of BH with M = 6 • 10""^^ g with 
respect to the total amount of DM in the Galactic Center 
is the same as the cosmological average, we would expect 
about 10^° Mq of DM in the Galactic Bulge. As we see 
later this is a reasonable result. On the other hand, we 
can avoid this restriction if the bulk of DM are not BHs 
but is of some other form, e.g. massive stable elementary 
particles. In this case the fraction of BHs in the Bulge can 
be larger than the cosmological average and the amount 
of DM in the Bulge may be significantly smaller than 
10^° Mq. 

Keeping an open mind to this possibility, we will con- 
centrate on a more economical case that all cosmological 
DM consists of primordial BHs, considering a more real- 
istic mass distribution than a delta function. The mass 
distribution of such BHs is model dependent and in par- 
ticular strongly depends upon the BH production mech- 
anism in the early Universe. One possibility is the model 
suggested in ref. [l^l, where in the simplest case the BH 
mass distribution has the log-normal form [s^] 



— = Cexpf-7ln ^ 



(10) 



where C, 7 and A/q are unknown parameters of the un- 
derlying theory. Let us note that this distribution is 
rather sharply peaked near mass A/q for 7 > 1. Strictly 
speaking the distribution given by eq. (jlOp is the ini- 
tial mass distribution, and the primordial BH masses 



changed in the course of cosmological evolution. There 
are two possible effects: evaporation which leads to mass 
loss and accretion which results in the mass increase. The 
latter leads to mass rise proportional to the square of 
the initial mass and weakly change the distribution. The 
evaporation results in complete destruction of primordial 
BHs with small initial masses, M < 5 • 10"'^^ g, and thus 
to decrease the spectrum at small masses. For the taken 
here mass distribution with Mq ~ 6 • 10"'^^ g and 7 > 1 the 
present mass distribution is quite close to the primordial 
one, because the mass fraction of BHs which have already 
evaporated is small (see below for more details). 

The total amount of DM in the Galactic Center is 
poorly known. The observed rotation curve requires that 
the total mass in the Galactic Bulge is roughly 10^° Mq 
and though this can be explained by the known baryonic 
components, see e.g. ref. [28|, it is natural to expect that 
the amount of DM is of the same order of magnitude as 
that of baryons. N-body simulations predict some cusp 
of DM in the center of galaxies [1^ , since the typical DM 
density profile is p ~ with = — 2. In conclusion, 
^ 10^" Mq is a reasonable upper bound on the amount 
of DM in the region where most positrons annihilate. So, 
if we assume for example that the total mass in the form 
of primordial BHs in this region is 5 • 10^ Af©, we find 
that the choice 7 = 1 and Afo = 6 • 10^^ g can explain 
the 511 keV line flux from the Galactic Center and be 
consistent with the observed galactic and extra-galactic 
gamma backgrounds. This can be seen as follows. By in- 
tegrating over all the masses, we flnd the normalization 
constant C 



dN 
dM' 



■dM 



7 = 1, J\/(,=6-10"5 1 



5 ■ 10^ Mp 







(11) 
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FIG. 2: Gamma ray spectra from primordial BHs (dark-blue 
solid curve) and of the measured background (light -red solid 
curve) from the Galactic Bulge in 7 cm~^ s~^ as a function of 
energy E in MeV. The BHs are assumed to have log-normal 
mass distribution, eq. (|10p . with the parameters: 7 = 1 and 
Mq = 6 ■ lO^*' g. The number of BHs is now normalized by 
the condition that their total mass in the innermost 0.6 kpc 
is 5 • lO'^ Mq. The gamma flux from primordial BHs does 
not exceed the ±25% uncertainty of the measured gamma 
ray flux (red dashed lines) and can produce enough positrons 
to explain the 511 keV line. 



4 




2 3 
E (MeV) 

FIG. 3: Cosmic isotropic gamma ray spectra from primor- 
dial BHs (dark-blue curve) and of the measured background 
(light-red curve) in 7 cm~^ s~^ sr~^ as a function of energy 
E in MeV. Here we assume that the BHs make the whole cos- 
mological dark matter and have log-normal mass distribution, 
10|) . with the parameters: 7 = 1 and A/q = 6 • 10^® g. 



eq 



and then we can obtain the photon flux arriving on the 
Earth 



1 



dN^ dN 



dE 47rr2 / dE dt dM 



dM . 



(12) 



Here we can use the initial mass distribution, eq. pop . 
and integrate out from to infinity because the correc- 
tions due to BH evaporation arc neghgible. dN^/dEdt 
is provided by eq. Fig. [3 shows the gamma spec- 

tra from the Galactic Bulge produced by the primordial 
BHs (dark-blue solid curve) and the observed gamma 
background of eq. ^ (red solid curve). The ±25% un- 
certainty of the gamma background (light-red dashed 
curves) is at least a reasonable estimate of the 2a curve, 
see refs. 

The isotropic cosmological gamma background pro- 
duced by primordial BHs in the whole history of the 
Universe depends only on the parameter 7 and Mq, 
if the normalization constant C is found from the re- 
quirement that the primordial BHs make the whole cos- 
mological DM, whose energy density is known to be 
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FIG. 4; Diffuse cosmic neutrino spectrum from primordial 
BHs in 1/ cm~^ s~^ sr~^ as a function of energy E in MeV. 
The assumptions are the same as of Fig. |31 



PDA! ~ 2.4 • 10^^" g cm^^. The estimate of the dif- 
fuse photon flux is (see e.g. ref. [30]] for the derivation of 
the formula) 



dE 



1 

An 



dz 



Ho h{z) 

(13) 

where dn/dM is the comoving number density of pri- 
mordial BHs with mass M, h{z) — y^Qmi^ + z)^ + 
and Hq the Hubble parameter today. The spectrum is 
presented in Fig. [3l together with the measured extra- 
galactic continuum [3lj 



dE 



B 



E 



1 MeV 



-2.38 



(14) 



where B = 6.4 ■ lO^"^ 7 cm"^ s^^ sr^-^ MeV^^ is the 
normalization factor. Eq. (|14p fits quite well the ob- 
served spectrum in the energy range 0.1 — 10 MeV. Con- 
trary to the simplest case of primordial BHs with equal 
masses, i.e. with the delta-function spectrum, primor- 
dial BHs with sufficiently wide spectrum may explain the 
observed annihilation line, be in agreement with the ob- 
served gamma background, and make all cosmological 
DM. 

Fig. m shows the expected cosmological neutrino back- 
ground from the evaporating primordial BHs. Unfortu- 
nately, such neutrinos are difficult to observe, because 
they are below the solar and terrestrial neutrino fluxes. 

One could also wonder about the impact of the pri- 
mordial BH positron and electron emission spectra on 
the diffuse Galactic backgrounds. Simple considerations 
suggest that MeV positrons and electrons produced in 
the Galactic Bulge are confined to the Galactic Bulge [llj 
and therefore their fiux on the Earth is negligible. Even 
if we take into account the fact that such primordial 
BHs could make the whole cosmological DM, and thus 
some of them are also in the neighborhood of the Solar 
System, one cannot expect any effect on the measured 
Galactic background, because most of the positrons are 
non-relativistic and are likely unable to reach a detector 
which is located on a satellite orbiting around the Earth. 

The positrons produced by primordial BHs have low 
energies and there is no contradiction with the bound de- 
duced in ref. Q . Actually, the constraint of 3 MeV is ob- 
tained assuming that all the positrons have the same en- 
ergy, so it is not easy to apply the result to more general 
cases. However, we can quickly see that in the scenario 
proposed in this paper the positrons with energies above 
3 MeV are about 5% of the total number of positrons pro- 
duced by Hawking evaporation. This is found calculating 
the integral 



dM M- 



,dN 
dM 



dE 



3 MeV 



Ey/E^ - ml 
1 -f exp(£;/T) 



(15) 



and dividing by the same integral evaluated between Wg 
and infinity. Including the grey-body effect, the correc- 
tion to this estimate is irrelevant, close to a factor 1.1. 
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FIG. 5: Gamma ray spectra from primordial BHs (dark-blue 
solid curve) and of the measured background (light-red solid 
curve) from the Galactic Bulge in 7 cm~^ s~^ as a function 
of energy E in MeV. The BHs are assumed to have all the 
same mass: M = 10^^ g (upper panel), M = 10^^ g (central 
panel) and M = lO^'^ g (lower panel). The number of BHs 
is normalized by the condition that they produce the right 
amount of positrons to explain the observed 511 keV line. 
Red dashed lines are the ±25% uncertainty of the measured 
diffuse gamma flux. 



The mass distribution of eq. (|T0|) includes also very 
light BHs which have already evaporated. For our choice 
7 = 1 and Mq = 6 • 10^^ g, the initial fraction of BHs 
with masses smaller than 5 • 10^"* g, and hence lifetimes 
shorter than the present age of the Universe, is negligi- 
ble, at the level of 10~^^. However, the total amount 
of matter converted into rclativistic particles from the 
birth of the Universe up to today is much larger: the 
dominant contribution comes from BHs which still exist 



in the Universe and are basically evaporating at a con- 
stant rate, because the process is slow and their mass 
is still very close to the initial one. Assuming that pri- 
mordial BHs make the all DM, today the cosmic energy 
density of the particles produced by BH evaporation is 
roughly 10~^ Pdm- The photon flux is the one reported 
in Fig. El 

Lastly, we would like to note that the BH masses are 
strongly constrained by the gamma ray continuum from 
the Galactic Center. Indeed, neglecting any possible cor- 
relation between primordial BHs and DM, we can con- 
sider three simple cases, where all the BHs have the same 
mass: M = 10^^ g, M = 10^^ g and M = 10^"^ g. 
In order to explain the observed 511 keV gamma flux, 
the total mass of primordial BHs in the Galactic Bulge 
should respectively be about 7 ■ 10"' M©, 7 ■ 10^ M© and 
9 • 10^^° Mq. The produced galactic photon background 
for these cases is presented in Fig.O For low BH masses, 
the BH temperature is high and too many high energy 
photons are produced. For high masses, the positron 
emission is exponentially suppressed by the Boltzman 
factor and to explain the observed 511 keV line we would 
need a large number of BHs, with the effect that the low 
energy gamma spectrum is too intense. Hence, we can 
conclude that, assuming that the positron flux is pro- 
duced by primordial BHs, the BH mass distribution has 
to be peaked around 10^ g and actually rather sharply 
peaked near 6 • 10^^ g. 

Conclusion — We have considered evaporating pri- 
mordial BHs, as a possible source of positrons to gener- 
ate the observed photon 511 keV line from the Galac- 
tic Bulge. The analysis of the accompanying continuous 
photon background produced, in particular, by the same 
evaporating BHs, allows to fix the mass of the evaporat- 
ing BHs near lO'^^ g. It is interesting that the necessary 
amount of BHs could be of the same order of magnitude 
as the amount of dark matter in the Galactic Bulge. This 
opens a possibility that such primordial BHs may form 
all cosmological dark matter. The background MeV pho- 
tons created by these primordial BHs can be registered 
in the near future, while the neutrino flux may be still 
beyond observation. The significance of this model would 
be difficult to overestimate, because these BHs would 
present a unique link connecting early universe and par- 
ticle physics. 

After this paper was completed, we became aware of 
the paper [3^ where a similar idea was explored. How- 
ever, the authors of this work assumed much lighter BHs 
which is in contradiction with our results. 
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Erratum 

We overestimated the total annihilation rate in the 
Galaxy almost by factor 2. Instead of the value 

3 • 10*3 s-i , (16) 

used in the paper, the correct one is 



Our error was based on the statement of ref. [6] of our 
paper where it is written that "the true annihilation 
rate is 3.6 times larger than would be deduced from the 
0.511 MeV flux alone". The corrected number of the 
annihilation rate results in shift down of the theoretical 
curves in figs. 1 and 5 by factor 1.8 - 1.9. We thank 
Pierre Jean for indicating to this error. 



1.6 • 10-*^ s-i 



(17) 



